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ABSTRACT 


A method  is  mathematically  derived  for  the 
measurement  of  the  thermal  conductivity  of  gases  which 
is  independent  of  raalatioi,  and  thus  particularly 
suitable  at  elevated  temperatures.  The  method  employe 
n rtnnvent tonal  hot  wire  oaII  with  its  wire  heated  by 
a sinusoidal  alternating  current  superimposed  upon  a 
direct  current.  The  operating  conditions  under  which 
the  radiation  introduces  no  appreciable  error  are 
indicated  as  well  as  the  inherent  limitations  of  the 
method.  Optimum  cell  design  and  operating  conditions 
are  suggested. 
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I.  Outline  of  Proposed  Method  for  Determining  Thermal 
Conductivities  of  Oases  at  High  Temperature 

The  method  proposed  to  measure  thermal  conductivity 
at  high  temperatures  employs  a horn  wire  ceii.  in  whicn  the 
wire  is  heated  by  a variable  cyclical  current.  A square 
wave  of  direct  current  or  any  alternating  current  whose 
wave  form  is  known  could  be  considered.  However,  a 
sinusoidal  alternating  current  superimposed  upon  a con- 
stant direct  current  has  been  assumed  in  the  following 
nathematical  analysis.  This  wave  pattern  was  chosen, 
not  only  for  mathematical  convenience,  but  also  because 
the  performance  of  electronic  equipment  is  best  known 
in  terms  of  sinusoidal  alternating  current. 

The  principal  reason  for  the  mathematical  analysis 
was  to  obtain  the  theoretical  calibration  equations  and 
to  determine  the  optimum  cell  design  and  operating  con- 
ditions such  that  the  response  is  sensitive  tc  varia- 
tions in  thermal  conductivity  in  the  range  expected  and 
insensitive  to  wire-to-wall  radiation.  However,  3lnce 
it  most  likely  will  prove  necessary  to  use  empirical 
calibration  to  account  for  unavoidable  errors,  any 
repetitive  current  ehape  could  in  practice  be  employed. 

The  method  developed  in  this  paper  is  based 
upon  the  fact  that  the  temperature  of  the  wire  will 
oscillate  above  and  below  its  average  temperature  cut 
of  phase  with  the  input  power  due  to  the  fact  that  the 
gas  has  a finite  thermal  conductivity,  specific  heat, 
and  density  and  thus  does  not  instantaneously  remove 
the  heat  generated  within  the  wire.  Therefore  any 
physical  or  electrical  property  which  is  a function  of 
this  phase  lag  may  in  theory  be  measured  to  determine 
the  thermal  conductivity  of  the  gas.  The  following 
mathematical  analysis  applies  this  principle  and 
indicates  the  properties  which  may  be  measured. 
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II,  Mathematical  Analysis  of  Differential  Equation  and 
Accompanying  Bounda cy  Conditions 

I 

A wire  of  radius  r-j_  mounted  coaxially  in  a cylindri- 
cal Jacket  of  radius  :r2  is  heated  with  superimposed  A„C. 

and.  D.C.  currents.  Heat  is  removed  from  the  wire  both  by 
conduction  through  the  gas  which  fills  the  Jacket  and  by 
radiation  from  the  wire  to  the  wall.  In  order  to  find 
the  temperature  of  the  wire  as  a function  of  time,  the 
following  assumptions  are  made: 

1.  The  wire  is  at  a uniform  temperature,  t^, 

at  any  instant  (true  for  small  wires  and  low  frequencies). 

2.  All  physical  properties  are  Independent  of 
temperature  (true  for  small  temperature  fluctuations). 

3.  The  external  Jacket  is  maintained  at  a constant 
uniform  temperature,  t2  (true  for  a heavy  Jacket  at 

steady  state). 

4.  Heat  radiated  by  the  wire  is  independent  of 
time  and  proportional  to  the  average  power  imput  to  the 
wire  (true  for  reasonably  small  net  radiation  and  tempera- 
ture fluctuations  small  compared  to  the  wire-to-Jacket 
temperature  difference.  It  is  to  satisfy  this  condition 
that  D.C.  is  required  as  well  as  A.C.) 

For  any  desired  test  conditions  the  extent  of  deviation 
from  each  of  these  assumptions  can  be  computed. 

The  differential  equation  that  applies  in  the  gas 
phase  is,  in  cylindrical  coordinates: 

= 1st  r,  (1) 

br^  r br  a b<3  1 d 

The  boundary  conditions  at  the  wire  and  Jacket  sur- 
faces are  obtained  by  heat  balances. 
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Thus  at  r = r. 


Rate  of  total  heat  int©  wire:  (l  + I0  sin<o©)  R 

Rate  of  heat  absorbed  by  wire:  ttt-^c^p'L 

, . bt 

Rate  of  heat  conducted  Into  gas:  -2TTr-,kL  "Sr 

? Tn‘ 

Rate  of  heat  radiated  to  Jacket:  yQ*  = y (lc  + — S- 


•)R 


Equating  the  heat  Input  to  heat  absorbed  plus  heat  trans- 
ferred yields:  2 

p 2 ^ Aq 

(l  + IQ  8in<0<S)T*  = TTr1  c‘p'L  V5  - 2TTr1kL  br  + y(1c+T“)R  (2) 

A similar  heat  balance  at  the  Jacket  wall  yields: 


- Y)(l"  + 


I- 


2 


6t 


)R  = - 2nrpkL  br 


(3) 


The  solution  of  the  differential  equation,  using 
the  previously  stated  assumptions  and  boundary  values 
is  presented  in  the  Appendix.  The  following  expression 
is  obtained  for  the  temperature  of  the  wire  as  a func- 
tion of  time: 


*1  - t2  + 


(l2+I2)(l-y)R 


2TTkL 


In 


rl  ur^^ic'p’L 


cos  (2^  + 0)  + 


I 


£4  I 

U i 


L' 


2p2  lIRE’j 

p 

co  c1  p'L 


cos  (w©  +e) 


(4) 


where  J2  I = ID  and  E^,  E^»  0, fc  are  dimensionless 

quantities  which  are  complex  functions  of  the  geometry 
of  the  apparatus  and  the  physical  properties  of  the  wire 
and  gas.  The  relationship  between  the  variables  can  be 
expressed  in  terms  of  the  following  or  equivalent  di- 
mensionless groups: 


It  is  important  to  note  that  none  of  these  dimen- 
sionless ratios  ie  a function  of  radiation,  -which  only 
enters  in  the  coefficient  of  the  log  term  in  equation  4. 
Therefore  the  measurement  of  any  one  of  these  four 
quantities,  together  with  the  knowledge  of  the  required 
physical  properties  and  the  dimensions  of  the  apparatus, 
should  yield  the  value  of  a of  the  test  gas  without  any 
corrections  for  wire  to  wall  radiation,  provided  the  cell 
is  operated  so  that  the  previously  stated  assumptions  are 
valid. 


III.  Measurement  of  the  Magnitudes  and  Phase  Angles 
of  the  Temperature  Waves 


An  electrical  method  may  be  derived  for  the  measure- 
ment of  the  above-mentioned  dimensionless  groups,  E-^, 

E2i  0y  € , in  the  following  manner.  In  the  previous 

derivation  it  was  assumed  that  the  resistance  of  the 
platinum  wire  used  in  the  cell  was  constant.  If  it  is 
now  assumed  that  the  wire  resistance  is  a linear  func- 
tion of  temperature  over  the  range  of  wire  temperature 
fluctuation,  the  resistance  of  the  wire  can  be  expressed 
aa  R = Rc(l  + st-^),  where  Rq  is  the  extrapolated  value  or 

R at  t = 0°C.  The  voltage  drop  across  the  wire  becomes 
E = (1  + I0  ein^9)R  = (i  + IQ  slnw9)R0(l  + at^  (9) 


and  combined  with  equation  (4-)  yields 


E = (iR0  + J2  IRo  sin  <*}©)  (l  + st2  + ln 


J2  sR  Rol^;  , 21I^rr  sEp  . 

d sin  (w  0 + 0)  - 0 . £ sin  %■» 


2 

2rrr]_«oc 1 p 


■nr^iic  * p 'L 


+ jg  sln(3ue  + 0)  + aiV;  gln(Sue  +6)  (1Q) 

Pttt^uC  * p 'L  ■nri^c,p,L 


- —xm amit  Jt *w- 


Therefore,  the  voltage  drop  across  the  wire  contains  a 
d.c.  component,  a fundamental  frequency  component,  and 
second  and  third  harmonics.  Consequently 
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- E2™>l2o’p’Z.«  /'I  + (n) 

2 “ilrP ) 

E"  = E3  J^r^o'P'Wi  + t\  (12) 

3 ^ j 


where  is  the  amplitude  of  the  second  harmonic  voltage 
component,  E -z  Is  the  amplitude  of  the  third  harmonic 

voltage  component,  and  6 and  0 are  the  phase  angles  of 
the  second  and  third  harmonic  voltage  components,  re- 
spectively, f Is  the  average  temperature  of  the  wire 
and  R Is  the  resistance  of  the  wire  at  t. 

Therefore,  In  principle,  the  measurement  of  the 
amplitude  of  the  second  or  third  harmonic  voltage  com- 
ponents, or  their  corresponding  phase  angles,  will  give 
values  of  thermal  conductivity  independent  of  wire-to- 
wall  radiation.  It  can  be  shown  that  the  phase  lag 
plus  the  phase  angle  (0  ore  ) equals  tt  radians. 

Initially  It  was  believed  permissible  to  use  a 
pure  sinusoidal  input  with  no  superimposed  direct 
current.  In  this  case  1=0  and  the  second  harmonic 
component  of  voltage  Is  absent.  For  this  reason  only 
the  equations  for  E'i  and  0 given  in  Appendix  C were 
evaluated  as  functions  of  the  dimensionless  groups 
listed  in  equations  (5)-(S).  The  values  were  computed 
to  ten  significant  places  using  the  IBM  Electronic  Data 
Processing  Machine-Type  7^1 • The  significant  portion 
of  the  results  is  tabulated  In  Tables  I and  II  and 
plotted  in  Figures  1 through  7»  (The  tables  and  the 
figures  are  found  In  the  Appendix.)  However,  it  may 
prove  preferable  to  employ  the  second  harmonic  due  to 
lesser  build- up,  In  the  required  electrical  circuits. 

In  this  case  an  entirely  empirical  calibration  could  be 
used. 


Figure  5 shows  that  E?  Is  most  sensitive  to  changes 

2 ^ 

In  a at  values  of  (Brp)  near  4 at  the  practical  ratio 

O 

of  r^/ro  equal  to  0.007.  At  values  of  (Br^)  above 
approximately  2 the  initial  assumptions  are  reasonably 
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valid.  It  is  expected  that  E^  will  he  less ,(|suited  for 
this  purpose,  although  no  calculations  of  E2  have  been 

made.  The  initial  A.C.  current  will  undoubtedly  have  a 
third  harmonic  component,  which  will  of  course  have  to 
be  filtered  out  of  the  A.C.  before  it  is  applied  to  the 
hot  wire.  This  third  harmonic  will  be  useful,  however,  as 
a standard  from  which  to  measure  changes  in  0 between 
calibration  and  test  runs. 

IV.  Analysis  of  Calculated  Data  and  Graphs 

The  dimensionless  group  E^  and  the  phase  lag  A = 
tt  - 0 were  calculated  at  various  values  of  pr2,  r^/r 2> 

and  c'p'/cp  (where  P = ii ) and  plotted  as  the 
logarithm  of  either  ^3)^  - E’^j  or  (Amax  - A) 

versus  Pr2  at  constant  values  of  r]y'r2  and  c’p'/cp. 

In  addition  advantage  was  taken  of  the  fact  that  both 
E^  and  A,  as  well  as  their  slopes,  are  zero  at  Pr2 

equal  to  zero.  The  curves  in  the  range  in  which  pr^ 

is  less  than  0.01  is  dotted  in  order  to  indicate  their 
approximate  shape,  since  no  values  of  the  required 
Bessel  Functions  were  available  in  this  region.  In  any 
event  it  would  not  be  desirable  to  use  such  low  values, 
because  it  would  mean  an  extremely  small  wire  radius 
and/or  a very  low  frequency,  and  consequently  a large 
wire  temperature  fluctuation  would  result. 

The  following  conclusions  have  bean  reached  by 
consideration  of  the  data  and  accompanying  graphs: 

1.  As  Pr^  increases,  E^  approaches  1/2  as  a 

limit  and  A approaches  the  value  or  tt/2  radians  (90°). 

This  is  also  the  case  as  c’p’/cp  approaches  cao  . 

\ 

2.  For  elevated  temperatures  and  normal  or-  low 
pressures  (high  values  cf  c’p’/cp),  E^  varies  sufficiently 

with  Pr2  to  be  suitable  for  measuring  k.  In  particular, 

this  response  may  be  used  for  gases  at  high  pressures 
and  moderate  temperatures. 

3.  It  is  expected  for  elevated  temperatures  that 
the  phase  lag  will  be  sufficiently  sensitive  in  the 
range  of  Pr2  between  1.5  and  3>  and  ri/r2  equal  to 
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approximately  0.007,  to  measure  the  thermal  conductivity 
within  one  or  two  per  cent.  This,  however,  would  mean 
that  for  most  common  gases  the  fundamental  frequency 
must  be  fairly  small.  It  has  been  estimated  that  60 
cycle  alternating  current  should  be  useful  for  air 
if  interference  from  power  lines  is  not  objectionable. 


V.  Conditions  for  Evaluating  Optimum  Ratio  of  Direct 
Current  to  Superimposed  Alternating  Current 

From  the  practical  viewpoint  of  experimental 
design  and  operation,  the  most  convenient  method  would 
be  the  limiting  case  of  a pure  sinusoidal  current 
with  no  superimposed  direct  current.  However,  it  has 
been  calculated  for  this  case  that  at  optimum  operat- 
ing conditions  the  temperature  fluctuation  of  the  wire 
is  about  30  to  oO  per  cent  of  the  average  temperature 
difference  between  the  wire  and  the  outer  Jacket. 

This  would  mean  that  the  assumption  that  instantaneous 
radiation  in  any  run  was  constant  and  proportional  to 
the  average  power  input  to  the  wire  wsb  not  valid  and the 
previously  mentioned  dimensionless  groups  would  have 
to  be  corrected  for  radiation  effects. 

In  order  to  decrease  the  temperature  fluctua- 
tion of  the  wire  to  a small  fraction  of  the  average 
temperature  difference  between  the  wire  and  outer  wall, 
direct  current  is  superimposed  on  the  alternating 
current.  The  restrictions  on  the  ratio  of  direct 
current  to  alternating  current  are  as  follows: 

1.  The  third  harmonic  component  of  voltage 

must  be  an  accurately  measureable  quantity.  Thus  there 
is  a minimum  alternating  current  input. 

2.  Because  it  is  assumed  that  all  physical 
properties  are  constant,  the  temperature  difference 
between  the  wire  and  outer  Jacket  must  be  as  small  as 
possible.  Therefore,  the  total  power  input  to  the 
wire  must  be  kept  small  by  limiting  the  amount  of 
direct  current  to  the  wire  after  the  A.C.  is  set  at 
its  operable  minimum. 

It  will  not  be  possible  to  satisfy  the  condi- 
tions of  a small  temperature  diff  eren  ce  and  a high 
ratio  of  direct  current  to  alternating  current  simul- 
taneously if  a minimum  value  of  alternating  current 


9 


Is  maintained.  Therefore  it  is  necessary  to  determine 
the  ratio  which  gives  the  most  reasonable  approach  to 
satisfying  these  restrictions.  It  is  expected  that  this 
optimum  ratio  of  direct  current  to  the  effective  alter- 
nating current  will  be  approximately  five  with  a 
maximum  average  temperature  different  of  about  40°C  and 
a minimum  effective  voltage  of  the  third  harmonic  in 
the  range  of  10“®  to  10“5  volts. 


V I.  Evaluation  of  the  Radiation  Assumption 


Due  to  the  Importance  of  the  assumption  that  radia- 
tion is  constant  in  any  run  and  proportional  to  the 
average  power  input  in  obtaining  the  result  that  E'-j, 

E*2»  €.  are  Independent  of  wire-to-wall  radiation, 

it  is  desirable  to  check  whether  this  assumption  is  valid. 
A more  vigorous  solution  of  this  boundary  value  problem 
would  be  obtained  by  recognizing  that  radiation  is 
temperature  dependent  according  to  the  Stefan-Boltzman 
law;  thus 


qr  = ^FAFE(2Ttr]L)(T14-Tp4)  = 


«■  FAFeV 


*1 
l + — 

T2 


(2-nr-jL)  (t-j-tg 


It  can  be  further  shown  that  if  the  temperature 
fluctuation  of  the  wire  is  2^  of  the  average  tempera- 
ture difference  between  the  wire  and  Jacket,  the 


assumption  that 


erW23 


[l  +f!i\  + f'l’NVVj5' 

L VTsi  v*i/  yj. 


is  a 


constant  and  equal  to  a heat  transfer  coefficient,  h, 
will  result  in  a maximum  error  in  the  value  of  qr  of 

about  0.15#.  The  heat  balance  on  the  wire  as  pre- 
viously stated 


Mfr  “N  - (i  + I0  sinw©)c  R - yQ»  (2) 
can  be  modified  for  this  more  rigorous  treatment  to  the 


6t 


) (13) 
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form 


M “N  |$  = (i  + IQ  sin  cj9)2R~yQ« 


tl~t2 

ti-t2 


(14) 


where (qp)=  yQ*  = hN(t  - t?)  and  = hN(t,  - tj 
avg  x r x 2 


If  equation  14  is  correct,  then  the  maximum 
error  in  the  left-hand  side  of  the  original  heat  balance, 
equation  2,  will  be  less  than  0.05#,  providing  that  the 
wire  temperature  fluctuation  is  about  2$>  and  that  the 
fraction  cf  the  total  power  input  transferred  by- 
radiation  (y)  is  less  than  1/5*  However,  because  this 
error  is  absorbed  by  every  term  involving  r and  0,  then 
it  is  expected  that  the  error  associated  with  either 
the  phase  lag*  A,  or  the  amplitude  of  third  harmonic, 

E^,  is  negligible. 


In  any  event,  if  it  does  become  necessary  to  make 
some  correction  for  radiation,  and  since  h and  the  true 
average  r-^  may  not  be  accurately  determinable,  an 
empirical  correction  for  the  given  apparatus  would 
probably  be  more  satisfactory  than  a theoretical  one. 
From  the  mathematical  analysis,  it  can  be  shown  that 
the  radiation  correction  will  be  a function  of 
hr-^/k  in  addition  to  the  previously  mentioned  dimen- 
sionless groups.  Substitui lag  in  the  definition  of  the 
heat  transfer  coefficient  in  terms  of  the  variables  com- 
bined in  the  Stefan-Boltzman  relationship  which  are 
equivalent  to  it,  yields  the  following  relationship 
for  the  radiation  correction: 


T I*! 

Radiation  Correction  = F i(3ro)  . 

2 °P  ' T~'  k 


(15) 


711.  Heat  Capacity  of  the  Gas 

Throughout  the  derivation  the  heat  capacity  of  the 
gas  was  used,  without  indicating  whether  it  was  the 
heat  capacity  at  constant  pressure  (cp)  or  at  constant 

volume  (cy).  It  has  not  been  definitely  determined 
which  term  is  the  correct  one  to  use.  Thus  it  is 
planned  to  ascertain  this  by  using  blank  runs  at  low 
temperature  with  gases  of  known  thermal  conductivity. 


However,  if  the  temperature  fluctuation  is  con- 
sidered to  move  in  the  gas  radially  in  the  form  of 
temperature  waves,  then  flow  waves  will  be  corres- 
pondingly moving  through  the  gas  in  the  same  manner. 
Since  these  fluctuations  are  occurlng  at  about  60 
times  per  second,  the  flow  waves  must  be  very  close 
together,  and  the  process  may  be  considered  to  be 
occurring  at  constant  pressure.  Therefore  it  is 
expected  that  Cp  will  be  the  proper  value  to  use* 


VIII.  Conclusions 


1.  It  is  shown  that  it  is  theoretically  possible 
to  obtain  the  thermal  conductivity  of  gases  at  high 
temperature  from  electrical  responses  which  are  suf- 
ficiently sensitive  to  thermal  conductivity  but  are 
insensitive  to  wire-to-wall  radiation.  A hot  wire 
cell  would  be  employed  in  which  the  wire  is  heated  by 
a sinusoidal  alternating  current  superimposed  upon  a 
direct  current. 

2.  The  phase  lag  of  the  third  harmonic  of 
voltage  produced  across  the  wire  with  respect  to  some 
standard,  such  as  the  third  harmonic  present  in  the  un- 
filtered generator  output, appears  to  be  a sensitive 
indicator  of  thermal  diffusivity  for  gases  at  high 
temperature. 

3.  It  also  seems  feasible  to  use  the  amplitude 
of  the  third  harmonic  to  measure  the  thermal  diffu- 
sivity of  gases. 

4-.  It  seems  possible  to  use  the  amplitude  of  the 
third  harmonic  to  measure  the  thermal  diffusivity  of 
liquids.  However,  calculations  at  the  low  volumetric 
heat  capacity  ratios  (c'p:/cp)  that  would  be  obtained 
with  liquids  have  not  yet  been  carried  out. 
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APPENDIX  B 


Solution  of  the  Differential  Equation 

The  differential  equation  and  the  boundary  condi- 
tions may  be  written  as  follows: 


it 

dr2 


dr 


W 


1 6t 
a d5 

where 

rl< 

r < r2 

(1A) 

- *yQ» 

= M $ 

- kN 

fct 

*r  - r - rl 

(2A) 

rl 

r„ 

at  r = 

r2 

where : 

(3A) 

M = nr^  c * p *L 
N = Pnr1L  and 

= ^ + 1o  sln 


2 /,  i02x| 

0 9)RaksMl  +“?“  Jr 


The  solution  of  this  boundary  value  problem  is  ob- 
tained by  assuming  a solution  In  the  form 

t = f + g cos  + h sin  2f-'9  + m sinuO  + %t  cosw©  (4-A) 

where  f,  g,  h,  1,  m are  unknown  functions  of  r and 
independent  of  6. 

Combining  equations  1A  through  4-A  yields  the  following 
set  of  equations 


f" 

+ 

1 

r 

f = 0 

(5A) 

% + f %'  = m 

(SA) 

s" 

1 

r 

g 1 = h£2 

(6a) 

m"  + i m*  = -Jfc 

(9A) 

il 

+ 

1 

r 

h*  = -gR2 

(7A) 

P2  = % 

(10A) 

Defining  the  following  sets  of  functions 

p = z ~ Jh  (llA)  q = % - Jm  (12A) 

where  j - (- 1 )3^2 

results  in  the  rewriting  of  the  differential  equation 
and  its  accompanying  boundary  conditions  as  follows: 


f"  + 

1 

f ' 

= 0 ' 

p"  + 

r 

1^ 

r 

P* 

= JP2p 

q"  + 

1 

r 

q* 

~ J T 

kNp  *- 

■2  J«Mp 

i02R 

i> 

JkNq*  +«*Mq  = -2iI0R 


r-,  <.  r < r. 


at  r = rn 


(5A) 

(13A) 

(lU) 

(15A) 

(16a) 


p»  = 0 
q ' = 0 


at  r = r2 


(17A) 

(IgA) 


Therefore  the  partial  differential  equation(lA) 
is  now  expressed  in  terms  of  three  ordinary  differential 
equations  (5A),  (13A),  and  (l4A),  which  are  easily 
solved,  yielding: 


-Q'CL-y)  . . _ 

f er  "X r~ In  v-  -f  C_ 

dTT/Jj  * O 

p = AJ^pJ^r)  + EK0(PJ 
r fl  3/2  \ 

q = J-  r;  + DK0 


A » \ 

1/2r) 

(20A) 

( & 1/?  \ 

(j?  J * 

(21A) 

CD,  A,  B,  C,  and  D are  constants  of  integration  depen- 
dent upon  boundary  conditions,  2A,  3A,  and  15A  through 
1&A*  The  constants  A,  B,  C,  and  D are  composed  of  real 
and  imaginary  parts,  Just  as  the  functions  p and  q ,by 
definition, are  complex  variables.  Thus,  to  separate 
these  functions  into  their  real  and  Imaginary  parts,  the 


I 
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Bessel  Functions  may  be  separated  in  the  following 
manner  ^7*  H) : 

3 *2  1/2 

Jo<2J  ' ) - her  z + J bei  zj  K0(zJ  )=  ^er  2 + 3 kei  2 

*x  / 2 “2^  1 /p 

J^(2J^'  ) = ber^  7.  v J bei^  zj  J K.^(zJ  )=kei^z+J  kej^z 

Accordingly  the  functions  p and  q can  be  separated  into 
their  real  and  Imaginary  parts,  as  indicated  in  the 
original  definition  of  these  functions  (Equations  11A 
and  12A)  and  may  be  re  substituted  into  the  original 
trial  function,  Equation  4-A. 

Finally  by  assuming  that  the  temperature  fluctua- 
tions  are  damped  out  at  r = and  therefore  that  is 

constant,  CG  will  equal  tp.  All  the  neceseary 

arbitrary  constants  are  thus  specified,  and  the  boundary 
value  problem  2?.  solved.,  An  additional  modification  is 
to  combine  similar  trignome trie  functions  into  one  func- 
tion by  the  use  of  a phase  angle. 
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APPENDIX  C 


Final  Result  of  the  Boundary  Value  Problem 

(i2+  I2)(R(1-y)  r2  I2RE^ 

o = tp  + Ztt'Aj  In  r-,  + 7~ cos(3w0+  0)  + 

^ x nr^'-e ' p 'L 


2 1 2 1ire_2  008  (we  + e ) 

ttt12c  * p 'L 

lo 

where  I = effective  current  = ~z^r 

12 


»•  2 2 

(E^)  « (A^ber  Pr-j_  •+  Ar,  bei  p^  - A-j  ker  pr-L  - A4  kel  pr^)  + 

2 

(A^  bel  pr-j  - A2  ber  pr-]_  - A^  kei  p^  + A^  ker  pr^) 

, * ,2  f pri  pr-,  pr-,  Pn 

(15.)  = lBl  ber  — - + B2  bei  -2  - B,  ker  — i - B4  kel  -=- 


r1  ™ vr 

- prj. 

B1  bei  ^=- B2  ber 


V* 

prx 

TT " 


is 

Ppi  + 


prx 


B-,  kel  y - 1 + Bk  ker 

J is  4 i* 


tan  0 « A^eipr-L  - AgberPri  “ A-lkelpr1  + A^kerp^ 
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APPENDIX  D 


Summary  of  Principle  Methods  of  Measurirg  Thermal 
Conductivity  of  Fluids 

The  following  principles  have  been  used  for  the 
measurement  of  thermal  conductivity  of  fluids: 

A.  Steady  State  Methods 

1.  An  electrically  heated  wire  in  a much  larger 
cooled  concentric  cylinder  filled  with  the  gas  or  liquid 
(2,  S). 

2.  An  electrically  heated  cylinder  in  a slightly 
larger  concentric  cylinder  filled  with  the  gas  or 
liquid  (5). 

3.  A disk  of  the  liquid  between  one  fluid-cooled 
plate  and  one  fluid-heated  (3)  or  electrically  heated  (l) 
plate . 

4.  Two  concentric  cylinders,  one  heated  and  one 
cooled,  of  which  at  least  one  is  non-radiant  and  non- 
absorbent (4). 

B . Unsteady  State  Methods 

1.  A constant  electrical  current  is  passed  through 
a thin  wire  in  a cylinder  of  the  liquid  initially  at 
uniform  temperature.  The  change  cf  temperature  with  time 
of  a point  close  to  the  wire  is  measured  (12). 

2.  A tube  containing  the  gas  is  electrically  heated 
by  a sudden  flow  of  electricity.  The  pressure  of  the 
gas  is  recorded  as  a function  of  time  (6). 

3.  A tube  ie  heated  by  an  alternating  current  of 
very  low  frequency,  and  the  time  is  measured  for  the 
temperature  wave  to  reach  the  center  of  the  viscous 
liquid  within  the  tube  ( 13 ) - 

The  steady  state  methods  are  mathematically  very  much 
simpler  than  the  unsteady  state  ones.  Methods  A-l  and 
A-2  are  the  most  popular  methods.  They  are  suitable  for 
gases  because  the  heat  losses  in  parallel  to  the  main 
conduction  through  the  gas  are  concentrated  at  the  ends, 
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can  be  minimized  by  using  long  cylinders,  and  can  be 
corrected  for  by  blank  runs.  Usually  they  are  not  con- 
sidered suitable  for  high  temperatures  because  of  the 
radiation  correction,  which  is  expected  to  be  high  as 
well  as  to  vary  with  time  due  to  change  in  surface  con- 
dition. Recently,  however,  a silver  cell  has  been 
employed  with,  a constant  radiation  correction  of  only  6% 
at  775 °C  (9).  Method  A-3  is  not  suitable  for  gases  be- 
cause of  the  relatively  greater  magnitude  of  the  losses. 

Method  A-4  constituted  an  initial  effort  to  develop 
a technique  for  high  temperature  gaseous  thermal  con- 
ductivity determination  which  is  free  of  the  wall-to-wall 
radiation  correction.  As  in  method  A-2,  two  concentric 
cylinders  arp  employed  with  the  test  gas  between.  How- 
ever the  inner  cylinder  is  constructed  of  a non-absorbent 
solid,  and  cooled  (or  heated)  by  a high  velocity  stream 
of  non-absorbent  gas.  Thus  no  net  flow  of  heat  by  wall-to- 
wall  radiation  exists.  In  actual  practice  a thin  fused 
quartz  inner  tube  was  employed  and  an  effective 
absorptivity  of  14J6  was  obtained  instead  of  a value  near 
zero  obtainable  with  the  more  fragile  KC1,  NaCl,  or 
GaF^  tubes.  Average  gas  temperatures  up  to  about  5°°°C 
were  reached.  At  higher  temperatures  the  effective 
absorptivity  would  have  been  even  lower.  End  effects 
were  determined  empirically  by  calibration  runs  with  air 
at  low  temperatures. 

The  unsteady  methods  are  in  general  all  quite  re- 
cent, having  been  published  subsequent  to  the  initiation 
of  the  present  study.  Method  B-l  is  not  suitable  for 
gases  because  the  thermocouple  would  influence  the  gas 
temperature  due  to  its  high  heat  capacity.  Method  B-2 
was  shown  (6)  not  to  be  reliable  because  of  its  de- 
pendence on  an  unknown  average  of  Cp  and  Cv.  Method 

B-3  is  only  suitable  for  highly  viscous  liquids,  having 
been  specifically  developed  for  molten  glass.  Further- 
more, the  values  obtained  are  uncorrected  for  radiation 
effects. 

Method  A-4  seems  to  be  the  most  suitable  for  the 
measurement  of  the  thermal  conductivity  of  gases  at  yet 
higher  temperatures  of  all  the  possible  methods  that 
have  been  employed  previously.  The  purpose  of  the  present 
study  is  to  investigate  another  possible  method  which 
will  measure  thermal  conductivity  without  interference 
by  wall-to-wall  radiation.  It  is  planned  to  push  forward 
the  application  of  whichever  method  seems  more  suitable. 


,a  Kc.  h-  p*  o ooo  ^ *-l  a Wltflioas^H 


NOMENCLATURE 
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I 


Voltage  drop  across  wire 

Third  harmonic  component  of  voltage 

Second  harmonic  component  of  voltage 

Heat  radiated  per  unit  surface  area  of  wire 

Geometric  Radiation  Factor 

Non-black  factor  of  surfaces 

Maximum  value  of  alternating  current 

Effective  value  of  alternating  current 
Length  of  wire  p 

Heat  capacity  of  the  wire  rrr^^c’p’L 
Surface  area  of  wire 
Average  power  Input  to  wire 
Resistance , electrical,  of  wire 

Resistance  of  wire  at  Its  average  temperature,  tn 

Resistance  of  wire  extrapolated  linearly  to  0°C  = 

R/fl  + at'1! 

* * 

Absolute  temperature  of  wire 
Absolute  temperate  of  wall 

Heat  capacity  of  gas 

Heat  capacity  of  wire 

Heat  capacity  at  constant  pressure 

Heat  capacity  at  constant  volume 

Heat  transfer  coefficient  for  radiation 
Direct  current  in  wire 

rr 

Thermal  conductivity  cf  gas 
Heat  conducted  through  gas 

Heat  radiated  from  wire  to  wall 
Radial  distance  from  center  of  wire 
Radius  of  wire 
Radius  of  outer  Jacket 

Change  In  resistance  with  temperature/RQ 
Temperature 

Temperature  of  wire  at  any  time 
Temperature  of  outer  Jacket 


t 


^•DT)  ii  '??,  (h  "CO  P cH 


2^4- 


Average  wire  temperature 
Thermal  diffusivlty  of  the  fluid 

l)F 

Fraction  of  power  Input  dissipated  by  radiation 
Phase  lag  of  third  harmonic  in  wire  voltage 
behind  the  fundamental  or  the  current 
Phase  angle  for  the  fundamental  of  the  tempera- 
ture fluctuation  with  respect  to  the  current 
Phase  angle  for  second  harmonic  of  temperature 
fluctuation 
Time 

Density  of  x'ne  fluid 
Density  of  wire 
Etofan-Eoltzman  Constant 
Angular  velocity  = 2tt»  (frequency) 


Dimensionless  Groups 


E-j  JPrra’i^c'p^ci 


E2TTri2c'p,L  w 
2iI2R2 
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